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The resistivity in metals near an antiferromagnetic quantum critical point (QCP) is strongly affected 
by small amounts of disorder. In a quasi-classical treatment, we show that an interplay of strongly 
anisotropic scattering due to spin fluctuations and isotropic impurity scattering leads to a large 
regime where the resistivity varies as T a , with an anomalous exponent, a, 1 < a < 1.5, depending 
on the amount of disorder. I argue that this mechanism explains in some detail the anomalous 
temperature dependence of the resistivity observed in CePd2Si2, CeNi2Ge2 and Celn3 near the 
QCP. 
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In the last five years an increasing number of heavy- 
fcrmion metals near an antiferromagnetic (AFM) quan- 
tum critical point was shown to display striking devia- 
tions from conventional Fermi- liquid behavior A 
few of the dirtier systems (e.g. Q|) appeared to show the 
T 3 / 2 behavior of the resistivity, p(T), as predicted by 
the Hcrtz-Millis spin-fluctuation theory for such a QCP 
|s|-[7|. In other cleaner single-crystal systems, such as 
CePd 2 Si 2 , CeNi 2 Ge 2 and Celn 3 §§], p(T) varies as T a 
with exponents a between 1.1 and 1.5. It has been argued 
H that such a behavior signals a fundamental breakdown 
of Fermi liquid theory. 

In this letter we propose a simple theoretical explana- 
tion for this anomalous behavior which covers both dirty 
and clean systems. Our main result is that the resis- 
tivity anomalies mentioned above can be attributed to 
the interplay between quantum-critical AFM spin fluc- 
tuations |5||| and impurity scattering in a conventional 
Fermi liquid. In fact, the resistivity can be described 
in semi-quantitative terms in the context of the simplest 
semi-classical Boltzmann equation approach. Here, we 
will have nothing to say about the striking linear behav- 
ior of piT) near the QCP in CeCu6- K Au K [Q in which 
the neutron scattering experiments suggest a more un- 
conventional behavior of the spin-fluctuations P] . 

The temperature dependence of p(T) can already 
be understood from the following qualitative argument: 
First we recall that scattering off AFM spin-fluctuations 
is most effective near "hot lines", i.e., points at the Fermi 
surface connected by the ordering wave vector Q (see 
Fig. [I]) where gaps would open up in the antiferromagnet- 
ically ordered (metallic) phase. As explained by Hlubina 
and Rice the strong scattering near these lines is 
short circuited by the "cold" regions on the Fermi surface, 
where the scattering rates are small. In clean systems the 
latter dominate the transport and resistivity acquires the 
usual Fermi liquid behavior, p(T) ~ T 2 at sufficiently low 
temperatures. Impurity scattering leads essentially to an 
averaging of the scattering rate over the Fermi surface re- 



ducing the effectiveness of the Hlubina-Rice mechanism 
and emphasizing the role of the "hot lines" . The above 
line of argumentation can be made more quantitative: 
Near the "hot lines" the scattering rate 1 /ts of the quasi- 
particles due to the quantum-critical spin-fluctuations is 
linear in temperature tm/ts ~ t — T/T where T is a 
characteristic energy scale and tm a typical scattering 
time. The width of these "hot lines" is given by y/i (see 
below) . In the "cold regions" , we expect Fermi liquid 
behavior with tm/t~s ~ t 2 . Weak disorder leads to an 
isotropic scattering rate, X/r e \ — x/tm, where £ is a small 
dimensionless number measuring the effectiveness of im- 
purity compared to magnetic scattering (x~ 1 rs kpl for 
spin-fluctuation scattering in the strong coupling "uni- 
tarity" limit; I is the elastic mean- free-path). The con- 
ductivity <7 is proportional to the average of over the 
Fermi surface (FS) with l/r k = l/r e i + l/rg(k): 
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The two terms describe the contributions from "hot" and 
"cold" regions, respectively. For x < t 2 < 1 cold regions 
short-circuit the hot ones and from the second term in 
(0) we obtain p oc i 2 |pjp|| ; while for t < x < 1, the hot 
spots dominate and the resistivity 1 j a is proportional to 
x + t 3 / 2 . At intermediate temperatures, x < t < i/x, 
we expect a crossover regime, in which we can define an 
effective resistivity exponent a, in terms of the logarith- 
mic difference of Ap(T) = p(T) — p(0) at the crossover 
temperatures T-f = Tx and T% = Tyfx, 



In Ap(T 2 c ) - lnAp(Tf) _ Inx - Inx 3 / 2 
In T% — In Tf In yfx — In x 
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In a very clean system this crude estimate implies a 
nearly linear crossover behavior in temperature as mea- 
sured in clean samples of CePd 2 Si 2 and CeNi 2 Ge 2 ||||. 

We now proceed with a more precise argument based 
on the semi-classical Boltzmann equation treatment of 
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electrons interacting with spin-fluctuations and impuri- 
ties. The former are described by a theory above the up- 
per critical dimension and, due to the ohmic damp- 
ing of the magnetic excitations in a metal, is character- 
ized by a dynamical exponent z = 2. As a result, the 
spin-fluctuation spectrum can be modeled by HQ 



XqH = X-qM 
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where qo and T are characteristic momentum and energy 
scales; and £ is the AFM correlation length which, at 
the QCP, diverges as l/£ 2 w eg 2 (T/T) 3 / 2 §,§. For the 
purposes of our numerical calculations we set c = 1 (c 
does not influence the low-temperature properties). 
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FIG. 1. Near the transition to an antiferromagnet with or- 
dering vector Q, the scattering on the Fermi-surface is en- 
hanced along "hot lines" connected by Q. The width of this 
region is given by Afc « qo^/T/T. The out-of-equilibrium 
distribution $e of the quasi particles for a current parallel to 
Q is shown as a function of the azimuthal angle 9 for temper- 
atures ranging from t — (qo/kp) 2 (T/T) 2 — 1 (left) down to 
t = 1CP 4 (right) and for both a clean system (lower curves) 
and a small amount of disorder, x — 0.01 (upper curves). 

Our starting point is the Boltzmann equation with 
a quasi-particle distribution function /k = — 
<f>k(<9/£/<9ek) linearized around the Fermi distribution 
with a collision term 



dt 
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Here gf mp and g 2 s are transition rates for isotropic impu- 
rity scattering and inelastic scattering from spin fluctua- 
tions, respectively, and n° is the Bose function. Eqn. (Q) 
tacitly assumes that the spin-fluctuations stay in equilib- 
rium, an approximation valid if the spin fluctuations can 
loose their momentum effectively by Umklapp or impu- 
rity scattering. Instead of solving the Boltzmann equa- 
tion directly, it is instructive to consider a mathemati- 
cally equivalent variational principle . Following Hlu- 
bina and Rice jio) , we define a functional p of <f>k 
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The physical resistivity is given by the minimum of /o[$k] 
regarded as a functional of $k- At this minimum <i>k 
is directly proportional to the true distribution function 
in an electric field applied in the direction of the unit 
vector n. In the above expression only integrals over 
the Fermi surface enter, as we have already integrated 
out the perpendicular components of k by using J dk — 
J de^ § dk/vk where Vk is the velocity of quasi particles 
at the Fermi surface. 

In the following we simplify the discussion by (i) con- 
sidering a spherical Fermi surface and (ii) limiting our- 
selves to transport parallel to the ordering wave vector 
Q = (0, 0, 2kp cos8h), in which case <&k = &e is only a 
function of the azimuthal angle 9. The equations 9 = 9h 
and 9 = tt—9h define the "hot lines" shown in Fig. |l|. The 
geometry @ and precise value of 9h are not very impor- 
tant as long as one stays away from Oh = ("2fcp" or- 
dering) or 9h — tt/2 (ferromagnetic ordering) where our 
approach breaks down fq|; in our numerical calculations 
we use 9h — 7r/6. As in |10| we approximate the second 
term in Eqn. (§) by 2ff|7 [y] with I[y] » 7r 2 /(y(3y + 2tt)) 
and y = {T/T) (l/(<?o£) 2 + (q ± Qf/ll) which is asymp- 
totically exact for large and small y. 

After performing the integration over u> and the polar 
angle tp in Eqn. (^) we obtain 

p(T) = rmn[p imp [<f> e }+ps[$e}} (7) 
xp M ffo ~ ®02 ) 2 sin M#i sin 8 2 d9 2 
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where pu — 3%|/(e 2 w|.) is a typical resistivity due to 
scattering from spin fluctuations at approximately the 
temperature scale T and x — .9 2 mp /(2<?f ) measures the 
relative strength of impurity scattering. The prefactors 
are chosen in such a way that xpM = 3fi.g 2 mp /(2e 2 t> fr) is 
the residual resistivity. The dimensionless scattering rate 
Fe 1 e 2 averaged over the polar angle tp is given by 
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$i = 9\ — 9h and d 2 = 8 2 — (it — 6h) measure the dis- 
tances from the "hot lines" and t = (T/r)(qo/kp) 2 is the 
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dimensionless temperature. Our numerical calculations 
use the full polar integral in ([l0|), even though at low 
temperatures, (l/(£kp) 2 ,t <C 1) the correct behavior is 
also contained in the approximate form. 

From Eqns. (0)-(^O|) one can easily deduce the qual- 
itative behavior of the low-temperature resistivity. We 
first consider the case of very low temperatures in a dirty 
metal. In this regime the resistivity is dominated by the 
disorder contribution p lmp . From dpi mp [Qe]/5Qg = 
we find the usual quasi-particle distribution for impu- 
rity scattering = cos 8. The leading temperature de- 
pendent correction to the residual resistivity po = xpu 
is then given by ps [cos 9] . The main contribution to 
ps arises in a small region around the hot spots, where 
- $e 2 ) 2 ~ (cosd H - cos(7r - Oh)) 2 is finite. Scal- 
ing $1/2 in Eqn. ( |l0| ) with \fi one recognizes that in the 
regime t > l/(kp0 2 ^ * the finite correlation length 
can be neglected and at lowest temperatures the resistiv- 
ity is given by 
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FIG. 2. Log-log plot of Ap — p(T) — p(0) for a rather clean 
system with x = 0.01. Note the large crossover regime from 
the resistivity of a clean system (dashed line) at high temper- 
atures to the resistivity of a dirty system (dot-dashed line). 
The inset shows, how this crossover evolves for various impu- 
rity concentrations x. 

On the other hand, if the system is clean (x = 0), 
we have to minimize pg . As pointed out by Hlubina 
and Rice |l0| , the ansatz <j># = cos is far from the true 
minimum. We can considerably reduce ps by using a 
distribution function where the "hot lines" are excluded, 
e.g. $ e = for \0 - H \ < cut and \9 - (tt - 9 H )\ < 
cu t- For this ansatz, the temperature dependence in the 
numerator of the scattering rate (|l0|) can be neglected 
for 9 cut ^> \ft and the resistivity is given by 



p(x = 0, T -» 0) = cp M (qo/k F ) 2 (T/rf 



(12) 



where c is a non-universal number of order 1 which de- 
pends on the details of the scattering mechanism in the 
"cold regions" of the Fermi surface. 



To obtain the crossover behavior we calculate the dis- 
tribution function and the resistivity p(T) within 
our model numerically by solving the integral equation 
5p[Qg]/6$g = which is equivalent to solving the lin- 
earized Boltzmann equation directly. For a clean system, 
<&e is shown in the lower part of Fig. [I]. At high tem- 
peratures the distribution function is structureless and 
all parts of the Fermi surface (besides those perpendic- 
ular to the current) contribute more or less equally to 
the resistivity. However, for lower temperatures the re- 
gion around the "hot lines" (dashed lines in Fig. |l|) are 
short circuited and the distribution function vanishes. 
Accordingly the resistivity is much lower and drops oc T 2 
(dashed line in Fig. ||). 
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FIG. 3. Effective exponent of the resistivity, defined as the 
logarithmic derivative of Ap(T). At very low temperatures, 
the "dirty-limit" exponent 3/2 is recovered. However, in the 
experimentally accessible low temperature regime smaller ex- 
ponents are to be expected for rather clean system (x < 0.1). 

As shown in the simple calculation discussed in the 
beginning, in a system with a small amount of disorder, 
we expect a large crossover regime between the behav- 
ior described by Eqns. (11) and (p^). In the variational 
approach given here this is due to the effect, that the 
impurity resistivity is not minimized by the distribution 
function $g lcan (the low temperature curve in the lower 
part of Fig. 0) and p imp [$f can ] = (l + c')xpM, where d is 
a number of order 1 (e.g. c' ~ 2.8 in our model). Below 
a temperature T|, defined by c'p(T = 0) « p{x = 0, T 2 c )i 
the distribution function deviates from $^ lean and ap- 
proaches the cos 0- form which minimizes impurity scat- 
tering (see Fig. |l|). Qualitatively, we obtain the same 
picture as in the crude estimate discussed at the begin- 
ning (up to factors like c'). 

The evolution of this crossover regime with impurity 
concentration is shown in Fig. ^| and its inset. The de- 
pendence of the distribution function (and of Ap) on 
impurity-concentration is a reflection of the complete 
breakdown of Matthiessen's rule in the crossover regime, 
where it is not possible to separate the different scatter- 
ing mechanisms contributing to the resistivity. In ad- 
dition, while not entirely physically meaningful, the IT- 
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dependent) effective exponent denned by the logarithmic 
derivative of p(T) — p(0) in Fig. |], when properly in- 
terpreted, displays the various crossovers in a dramatic 
way. For example, even for reasonably clean system ei- 
ther asymptotic exponents, 2 and 1.5 are difficult to ob- 
serve, while effective exponents close to 1 dominate over 
a wide range of parameter values as suggested by our 
estimate (@). 



40 



G 

a 
20 




0.08 



, 0.06 i 




0.00 



0.01 0.02 

:T/r(q /k F ) 2 ]' 



FIG. 4. Resistivity as a function of T ' in CePd2Si2 taken 
from g (left fi gure) compared to our calculation (right fig- 
ure) for x — 0.01. The insets show the corresponding loga- 
rithmic derivative of p(T) — p(0). The solid line in the inset 
of the theoretical plot displays the logarithmic derivative of 
p(T) -0.995p(0). Below « 400mK CePd 2 Si 2 is superconduct- 
ing (lower inset). Note the offset of the line T = in both 
plots. 

While the behavior for l/(kp^) 2 , t>\ are highly non- 
universal and strongly affected by details of the band 
structure and scattering mechanisms, this is not the case 
in the opposite limit l/(kp^) 2 , (< 1, which is of interest 
here. In the latter regime crossover effects depend only 
very weakly on the precise details of the model. 

We argue that our approach explains the anomalous 
resistivity at the QCP observed in the very good sam- 
ples of CePd 2 Si 2 , CeNi 2 Ge 2 and Celn 3 [§§]. In Fig. | 
we show the measurement by Grosche et al M of the 
resistivity as a function of T 12 at the critical pressure in 
CePd 2 Si 2 and compare them to our model for x = 0.01. 
It is important to note, that not only the same exponent 
of the resistivity shows up in the theory, but that it is 
also observed over a similar range O.lpo < Ap(T) < 10po- 
This suggests that standard spin-fluctuation theory 
can be applied for this system opening the possibility for 
a theory of the striking superconducting phase observed 
below 400mK. In CeNi 2 Ge 2 , a similar exponent is ob- 
served in the resistivity || , while the effective exponents 
in samples of Celn3 || show a behavior reminiscent of 
our predictions for x = 0.1 in Fig. |^. It is certainly nec- 
essary to check other predictions of this theory carefully. 
For example, according to ||, the pressure dependence 
of the Neel temperature near the QCP should be given 
by (p — ]3 C ) 2 / 3 while experimentally a linear dependence 
seems to be observed over some range || . Also the spe- 
cific heat should give valuable informations on the nature 



of the spin fluctuations. The most direct test of the ef- 
fects described in this letter is, however, a comparison 
of the critical resistivity in samples of different quality. 
According to our theory, the effective exponent has to 
change from 1.5 for dirty samples to values near 1 for 
very clean samples. Also, for cleaner and cleaner sam- 
ples, a "bump" has to show up in plots of the effective 
exponent (cf. Fig. ||). The dependence of the exponent 
on sample quality has indeed been reported 0. 

For systems not directly at but still near the QCP 
(kF^(T — 0) ^> 1), we expect again a large crossover 
regime in the resistivity with anomalous effective expo- 
nents due to a pronounced crossover from p ss pm(T/T) 2 
at high temperatures to p « pu{x + (T/P) 2 (£;f£)) at 
lowest temperatures. 
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